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,<1> : Abstract 

Q ■ We have calculated the leading order amplitude of hard exclusive production of two 

CN . pions in lepton nucleon scattering. At leading twist a pion pair can be produced only 

^ . in an isospin one or zero state. We have shown that isoscalar states are produced 

^ ! predominantly for xbj > 0.3 and with an invariant mass of the two pions close to 

the threshold (S-wave) and in the /2 resonance region (D-wave). These isoscalar 



O 



cn ■ pion pairs are mostly produced by two collinear gluons. Hence, comparing the 

^ ■ production of charged and neutral pion pairs as a function of xbj and m^^r, one can 

Q\ . get information about the gluonic component of two-pion distribution amplitudes. 
0^ 



Introduction 



{iLii Hard exclusive reactions open a new way to study the partonic structure of hadrons. The 

^ ■ factorization theorem []I| states that the amplitude of the reaction: 

-1, -fi + T^F + r (1) 

X. 

H , at large collision energy W ^ oo, large virtuality of the photon Q^ — »• cxo, fixed Q'^/W'^ 

(Bjorken limit), and with the masses M^, Mt', Mp <^ Q can be written in the form: 

Y^ j dz j dxifl/'j,{xi,xi - x;t, fi) Hij{Q^xi/x,Q'^, z, fi) ^f{z, ^) 
+power-suppressed corrections, (2) 

where fjij, is a T ^ T' skewed parton distribution |jl], |^, ^, ^ (for a review see 0), 
$J(z,/i) is the distribution amplitude of the hadronic state F (not necessarily a one 
particle state), and Hij is a hard part computable in pQCD as a series in asi^Q"^)- For 
the experimentally accessible range of Q^ and W the significance of this factorization 
might, however, be limited by the size of higher twist effects ||^, |^. Let us also note that 
often NLO corrections in hard exclusive reactions are noticeable |^, |^, ^. Therefore, in 
principle, reliable NLO calculations of the hard part Hij are desirable. 

In the present paper we shall study the hard exclusive electroproduction of two pions to 
leading order in the strong coupling constant as- The related process of photoproduction 
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Figure 1: Typical leading order diagrams for hard exclusive pion pair production. 



was studied in Ref. |11]. At leading twist the pion pair can be produced in a state with 



isospin one or zero. As argued in |T^ at small xbj, i-e. xbj < 0.1, the pions are produced 



mostly in the isovector state because the dominant reaction mechanism is exchange of 
two gluons with positive C-parity. Our prime interest will be the production of pions 
in the valence region of Bjorken x, i.e. xbj > 0.3. In this region (relevant for CEBAF, 
HERMES, COMPASS) the production of two pions is dominated by qq exchange, which 
leads to a sizeable admixture of pion pairs with isospin zero [T3\. We shall compute in the 
present paper the ratio of the cross sections for / = and 1=1 pion pairs as a function 
of Bjorken x and the invariant mass m.^T^ of the produced pions. 

The interesting feature of / = production of two pions is that these pions can 
originate not only from a coUinear qq pair but also from a pair of collinear gluons, see 
Fig. ^). Therefore measurements of J = hard pion production could give important 
information about the gluon content of the two-pion distribution amplitude. We predict 
the most favorable kinematic range for such measurements. 

Amplitude of the hard two-pion production 

In this section we compute the leading amplitude in l/Q"^ of hard two-pion electropro- 
duction. Due to the factorization theorem for exclusive hard reactions |l[] this amplitude 
can be written as a convolution of the skewed parton distributions in the nucleon (SPD), 
the distribution amplitudes of the produced pions (27rDA) , and a hard part computable 
in perturbative QCD as a power series in a.^- Hereby we consider the processes: 



7l{q) + B,ip)^27riq') + B2ip + A) 



(3) 



where A = p' — p is a four- momentum transfer to a target baryon Bi, —q^ = Q^ goes 
to infinity, and xbj = Q'^/'^P ■ Q is fixed. The state B2 can be a nucleon or a heavier 
baryonic state. Let us note that the corresponding amplitude for a transversally polarized 
virtual photon is 1/Q power suppressed. Therefore the considered (sub-) process with the 
longitudinally polarized photon, for which the factorization theorem holds, is sufficient to 
calculate the leading contribution of the whole electroproduction process e~{k) + Bi{p) —>■ 
27r(g') + e-{k-q) + B2{p + A). 

The leading order amplitude corresponding to diagrams of the type as shown in Fig. |T] 



has the formQ: 

{B,{p'),7rH\q')\r"^--SL\B^{p)) = -(e4na,)^j;^— j^^dr j^ dz 



X 



■f,f' ^ 



iV2 4Q 7-1 



^(r + O-^O (l-^)(r-0+^OJ 






(4) 



Here Cf is the charge of a quark of flavor f = u,d,s in units of the proton charge (e„ = 
2/3, Cd = Cs = —1/3), ^ is the skewedness parameter (see definition below), and ( is 
the longitudinal momentum fraction carried by the pion with flavor a. The longitudinal 
photon polarization vector is chosen as: 

£i = ^(g',0,0,g°). (5) 

The function Fffi (r, C,) is a skewed quark distribution defined as: 

r d\ 

FfAr,0 = J -e^'^{B,{p')\T{M->^n/2)hMXn/2)}\B,{p)) , (6) 

with n = 'Jan"' and the light cone vector n normalized so that: 



,2 



n 



0, n-{p + p') = 2. (7) 



*/'/(^'0 = / T7-e-^'^(^'"*\7rV|T{^;(An*)n>;,(0)}|0) 



The quark two-pion distribution amplitude ^J,j{zX) is defined as 

dX 

2^ 

and the two-pion gluon distribution amplitude ^q{z,() as: 

$«j'(z,C,mL) = -^/^e-^^('^'-*)(7rV|K'^n*'^G^^(An*)Gf (0)}|0), (9) 

where for convenience we introduced another light cone vector n* such that n ■ n* = 1 
and n* ■ n* = 0. The variable z is the fraction of the longitudinal (along the vector n) 
meson pair momentum q' = q — A carried by one of the quarks or gluons respectively. 
The variable ( characterizes the distribution of the longitudinal component of q' between 
the two pions: 

C=<^;-^ (10) 



^We write the most general result including the contribution of graphs of the type Fig. [|c) which are 
proportional to the skewed gluon distribution in the nucleon. It contributes only to the production of 
pions with isospin one and is very small in the valence region. Hence we shall neglect this contribution 
in the numerical estimates. For the skewed gluon distribution we use Fc{t, ^) in the notation of Ref. B. 



where k"' is the moinentuiii of the pion vr". 

The vectors n and n* are hnear combinations of the photon, initial, and final baryon 
momenta: 

M' = liMl+Ml-^). (11) 

The transverse plane is a plane orthogonal to the plane defined by the light-cone vectors 
n and n*. The skewedness parameter ^ (the longitudinal component of the momentum 
transfer) is defined as: 

^ = -lin.A). (12) 

In the Bjorken limit Q^ -^ oo, the skewedness parameter ^ is expressed in terms of Bjorken 

X as 4 = 2^- 

The two-pion distribution amplitudes (^J^) can be flavor decomposed as: 
TT+TT^ production: 

$i'4„(z,C,m^^) = (5^'^<l>^=°(z,C,m^0 + r/'V=i(z,C,m^0, (13) 

$^^""(z,C,m,,) = <l>^(;2,C,m,0, (14) 

^0^0 production: 

<^CoUz,C,m^^) = 6f'^<l>'=%z,C,m^^), (15) 

$^°-°(z,C,m^^) = <l>«(^,C,m^^), (16) 

where <l>^ are the leading twist quark 27rDAs with isospin /. For general properties of the 
two-pion distribution amplitudes see the next section. 

The leading twist skewed quark distribution with f = f = q = u,d,s and Bi = B2 = 
proton (Eq. (^) can be decomposed into nucleon spin non-flip and spin flip parts. Here 
we adopted the notations of Ji [Q for the spin decomposition of the matrix element of a 
bifocal quark operator between proton states. The spin non-flip part is denoted by Hf, 
the spin flip part by Ef. They are defined by: 

fix - 

—e''^{p'\^f{-Xn/2)n^f{Xn/2)\p) = HfiT,^,t) U{p') n U{p) 

+ ^ EfiT,^,t) U{p') a,,n^A^U{p) . (17) 

The independent variables of the skewed quark distributions Hg{T, ^, A^) and Eg{T, ^, A^) 
are chosen to be r, related to the fraction of the target momentum carried by the inter- 
acting parton Xi by: 

2 xi — X 
^" 2-x ' 

4 



the square of the four-niomentuin transfer A^ = t, and the skewedness parameter ^ defined 
in Eq. ([T^)^. In the forward case p = p' both A and ^ are zero and the second term on 
the r.h.s. of Eq. (|I^) disappears. The function H becomes the usual parton distribution 
function: 



if,(r,e = 0,t = 0) = g(r) 



-q{-r) 



r > 
r < 



Using the symmetry properties for the 27rDAs 
production can be written as: 



the leading order amplitude for tt+tt" 



{p{p').'n^Ti-{q')\r'^--eLW)) = -Mvra, 



N^. 



1 1 



N^ 24Q 

zil — z) 



U{p')nU{p) f^ rfr (2 //„(r, i, t) - H,{t, i, t) 



a r 



M a^.n^A' 



-f/(p)|_'^rfr(2E„(r,e,t)-i?d(r,e,t))a-(r)| 






+ 







zil-z) 



1 



+ U{V') 'l^^ U{p)j^^dr{2 E^{T,^,t) + E,{T,^,t))a+{r) 



(19) 



where 



a^(r) 



± 



T + ^-iO r-^ + iO 



(20) 



Let us note that for the production of two pions in an isovector state there is an additional 
contribution proportional to the skewed gluon distribution in the nucleon, see Fig. |I]c). 
The corresponding expression is given by the last term of Eq. (^ and also can be found 
in Ijl^, |T5[. Such a contribution is absent if the pions are produced in the isoscalar state. 
In the present paper we are interested in two pion production in the valence region and 
in particular in the isospin zero state in order to identify the gluons of Fig. ^). There 



the contribution from gluons in the nucleon can be safely neglected [jT6|. 
For the vr^vr'^ production amplitude we have: 



{p{p'),rr%'{q')\r-'^- ■ eL\p{p)) = -(e47ra. 



N? -I 1 



N^ 24Q 



X / dz 
Jo 



1 , (l-2z)<l>^=0(^,C,m,,; 



-^'^^i^X,m, 



zl-z) 



The skewed parton distributions and meson distribution amplitudes are scale dependent. The scale 
is set by the photon virtuality Q^. We do not show the scale dependence of these quantities to simplify 
the notation. 



X lu{p')hU{p) f^ dr{2 H^{t, e, t) - Hair, ^, t))a-{r) 

+ U{p' f'^^^ U{p)J_^dr{2 E^{r,i,t) - E,{T,^,t))a-{r)\ . (21) 

The expression for the ratio of the 1 = and 1 = 1 amphtude was obtained recently 
by Diehl et al. |TB[, however, the authors missed the contribution of the gluon 27rDA 
(Fig. [l|b)). Below we shall see that the contribution of this amplitude to the total two- 
pion amplitude is larger than the contribution of the quark 27rDAs. 

Two pion distribution amplitudes 

The two-pion distribution amplitudes describe the fragmentation of a pair of collinear 
partons (quarks or gluons) into the final pion pair [^ . As we shall discuss in the following, 
they can be related by crossing to the ordinary parton distributions in a pion. From C- 
parity one can easily derive the following symmetry properties of the 27rDAs: 

$^=0(z,C,m^^) = -$^=°(l-^,C,m^^) = $^=0(z,l-C,m^^), 
^^{z,C,m^^) = <l>^(l-2;,C,m^^) = $^(2;,l-C,m,,). (22) 



Following pD|, [T^ we decompose both quark and gluon two-pion distribution ampli- 



tudes in conformal and partial waves. For quark 27rDAs the decomposition reads: 

OO Jl+1 



$^=°(^,C,mL) = 6zil-z)j:T.Bifim^.)C'J'i2z~l)P,{2C-l), 

n=\ «=0 
odd even 

OO n+1 

^'=\zX.ml:) = 6z(l-z)EEi?r(^..)C^^/'(2^-l)PK2C-l), (23) 

n=0 l=\ 
even odd 

while for gluon 27rDAs we have: 

OO n+2 

$«(;^, C, mL) = 30 z\\ - zf Y. H Al{m^^)C'>J\2z - \)P,{2C - 1) . (24) 

even even 

Here Pi{x) are the Legendre polynomials and C^{x) are Gegenbauer polynomials. 

The distribution amplitudes are constrained by soft pion theorems and crossing rela- 
tions. The soft pion theorems for 27rDAs read ||T0| , p!2| : 



$^=°(z,C = l,m^^ = 0) = $^=0(^,C = 0,m^^ = 0) = 0, 
$'=^(z,C = l,m^^ = 0) = -<l>'=i(z,C = 0,m^^ = 0) = y^(2), 
<l'^(z,C = l,m^^ = 0) = <l>^(;2, C = 0, m^^ = 0) = . (25) 

Here ^{z) is the pion distribution amplitude. 

Additional constraints for the quark and gluon 2'7rDAs are provided by the crossing 
relations between two-pion distribution amplitudes and the known parton distributions 

6 



in the pion. For the derivation of such relations for quark 27rDAs see [O, for the gluon 



27rDA see [|^. The crossing relations have the form: 



^Af-l,Af(0) = ^ 



2 2N +1 



3 A^ 



1 1 

dx x^^'^—- ^{qn{x) + q^{x)) 



N. 



f q 



_ 2 2N + 1 /■! 

-^^'^\ Af(0) = o ^^7 — r / ^^ x^~^{u^{x) - ■u^(x)) , 



''Af- 



^N-2,nW 



3 N+1 Jo 



2N + 1 



(26) 



dx X g-wi^x) , 



5 {N + l){N + 2)Jo 

where Qnix), qnix), and gn{x) are the usual quark, antiquark, and gluon distributions in 
the pion. Using these relations one can easily derive the normalization of the two-pion 
distribution amplitude (27rDA) at m^^ = 0: 



0) = -— M^a^ - 



[ dz{2z-l)^^=^{z,C,m^ 
Jo 

Jo 

/' dz ^""{z, C, m^^ = 0) = -2 Mf C(l - C) • 
Jo 



(27) 



Here M2 and M^ are momentum fractions carried by quarks or gluons respectively in 
the pion and ^7^(771^7^) is the pion e.m. form factor. The soft pion theorems and crossing 
relations constrain $-'^"°(z, C,m^^) and ^'-^ {z , ( , ^TT-nn) only at 771,^^ = 0. For higher m^^ 
one can apply the dispersion relation analysis of Ref. [Iir 



The asymptotic expression for the isovector 27rDA has the form |T8 

<^i7y{z, C, m^n) =Qz{l- z) (2C - l)F^{m^^) . (28) 

Under evolution <l>*^ and <l>^=° mix with each other |T^. For asymptotically large Q^ one 

240 Cf 



obtains p!0[] 



^ ^^r,\Z^ C, , lll'J\ 



■ asy\ 



asy " 



z,C,m^ 



0)^ 
= 0) 



Nf + ACf 



120 



Nf + ACf 



z\i-zYa^-o, 

z{l-z){2z-l)C{l-C). 



(29) 



In the following we will use the asymptotic expressions for the 27rDAs. 
Numerical results and discussion 

Using our results for the production amplitudes discussed in the previous sections we can 
easily estimate the ratio of the cross sections for hard exclusive two pions production in 
the isoscalar and isovector states. The amplitudes of physically accessible reactions can 
be expressed in terms of / = and 1=1 amplitudes as follows: 






M^=° + M 



1=1 



(30) 



For the ratio 



R{xBj,m^TT,t) 



\M 



1=0 12 



\M^- 



:/|2 



(31) 



(integrated over the scattering angle of the pions) the xbj , t, and m^,r dependence factor- 
izes in the leading order: 



R{xBj,rnnn,t) = RDA{mTnr)RsPD{xBj,t) 

where the rrinn dependence is governed by the 27rDAs: 



RDAi^mr) 



with 



D^=^{m^ 



D^=Um^ 



(3 / dcosO 



D^=Hm^ 



(32) 



(33) 



P I dcos9 



1 _ (l-2z)$^="(^,C,m,,)-^<l>^(z,C,m, 

dz —^ 

z{l — z) 

1 $^=i(z,C,m, 
dz- 



zil-z) 



(34) 



Here 9 is the scattering angle of the pions in their cm. frame relative to the direction of 



g' 



/5cos^ = 2C-l with (3 



\ 



Ami 



mt 



(35) 



The dependence on Bjorken x = 2^/(1 + C.) and the squared momentum transfer t is 
determined by integrals over the skewed parton distributions (SPD): 



RsPD{xBj,t) 



\2i-{u)-iAUW 



|2/+(e,t) + /;(e,t)|2 



with 



lf{U) = j\dTHf{r,^,t) 



1 



± 



T + ^-iO r-^ + zO 



(36) 



(37) 



For the quark skewed parton distributions in the integrals Ir {^, t) we used the mod- 
els resulting from Radyushkin's double distributions [Q. They are obtained assuming 
F{x,y) = -K{x,y) ■ q{x) with the profile function 7c{x,y) = (^2^)3 suggested in pO[ and 



the MRS (A') parameterizations for the usual (forward) parton distributions q{x) ||21|| . 
The SPDs Hf(T,^) are related to Radyushkin's nonforward parton distributions JF/(X) 
by: 



Hfir,0 



1 + e 



0(r + O-^c 



e-r 



c'(fi|)-e«-)^/.i,, 



(3J 
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Figure 2: The shape of two-pion mass rriTrn (GeV) distributions for pions with isospin one 
(dashed curve) and isospin zero (solid curves). The isospin zero distributions are plotted 
for Bjorken x = 0.3,0.4,0.5. (The larger xbj the more enhanced is the distribution.) 



with ( = Y77, while the nonforward parton distributions are obtained from the double 
distributions by the integral 

j-c(x) = 



min{f,^} 



F{X-Cy,y)dy. 



(39) 



For the t-dependence we use the factorized ansatz H{T,C,,t) = H{T,^)Fi(t) with the 
proton form factor Fi{t). For simplicity we do not consider the specific contributions to 
quark SPDs which are not taken into account by double distributions. These contributions 



were observed first in chiral quark soliton model calculations |22| but actually they are 



of quite general nature as was shown in p3[. Also we neglected the contribution of the 
nucleon spin flip skewed parton distributions E. The complete studies, involving all these 
refinements, will be published elsewhere. 



Using the asymptotic expressions for the 27rDAs ( p8| , p9[ ) we can compute the integrals 
over z entering Eq. 



^.ii 



2z)^'=^{zX,m^^) 



40 



zil - z) 



Cp Jo z{l — z) 



Nf + ACf 



I T7^ h[m-n 



12 



Pn(cos( 



/3^ 



— -r f2{m^7r) P2(cos6')|, 



80 r3C-(3^ ^. 

4c;t^^ /oK 



1 ^^-\zX,m^^) 

dz ^ 

z{l — Z) 



Nf + ACf ^ 12 

77- /2("^7r7r) P2(cOs6')|, 

Q(3FJm^^) PAcosO). 



Po(cos( 



(40) 

(41) 
(42) 



(The effects due to deviations from tlie asymptotic forms were considered in |T2| , p^). We 
see tliat for tlie production of the pion pair in the isoscalar state the contribution of the 
gluon 27rDA is two times larger than the contribution of the quark 27rDAs. Therefore the 
hard production of the pions in the isoscalar state mainly probes the gluon structure of 
the pion. 

For our estimates we use the parameterization of the pion electromagnetic form factor 
F7r(m,r7r) glveu lu Ref. [^]. This parameterization reproduces the results of two- loop 
chiral perturbation theory at small 771,^^ and is in very good agreement with experimental 
data pH . 

The functions fo,2(}V), the so-called Omnes functions ^^, can be related to the vrvr 
phase shifts (5o(Vr) and ^^(Vr) using Watson's theorem |2^ and the dispersion relations 
derived in 



fiim^ 



exp 



z6? 



,m^ 



+ 



mz 



-Re 



TT 



ds- 



S?{s) 



4ml 



Sis — mt 



-iO) 



(43) 



The Omnes functions were analyzed in detail in Ref. |2^ and for the estimates of the 
cross sections one can use the results of that work. 

The constant C in Eqs. (|40|) plays the role of an integration constant in the Omnes 
solution to the corresponding dispersion relation. From the soft pion theorem it follows 
that C = 1 + 0{mD. Using the instanton model for calculations of Bniirn-^T,) at low 



energies [12, |3] one finds |10| the constant C to be equal to: 



C = 1 + 6 m^ + O(m^) with 6 ^ -1.7 GeV" 



The Omnes functions contain information about the tttt resonances as well as about 
the nonresonant background. Due to the /2(1270) resonance the /2(m,r7r) has a peak 
at mTTTT = 1.275 GeV. For the numerical calculation we use the Fade approximation for 
the Omnes function fo{Tn.„T,) suggested in Ref. fS^. For estimates of f2{jnT,T,) we used 
the expression ( ^31 ) where the D-wave pion scattering phase shift we saturated by the 
contribution of the /2(1270) resonance. 

In Fig. we plot the two-pion mass distributions in the isoscalar and isovector channel 
at different values of Bjorken x. 

The two-pion mass distribution in the isovector channel at large Q^ is determined by 



the pion e.m. form factor |T2 



dN 



1=1 



dm^ 



oc 36 f3 m^^ |F^(m7r7r; 



(44) 



which is plotted in Fig. ^ as a dashed curve. For the two-pion mass distribution in the 
isoscalar channel (for t = tmin) we can use: 



dN 



1=0 



dm-n 



J^y-^Bj) ^TT'TTTTy ^ri 



dN 



1=1 



drar, 



(45) 



dN'-° 



where R{xBj,'rn,,T,,t) is given by Eq. (^T]). In Fig. we plot ^£^ — (solid lines) for xbj = 
0.3, 0.4, 0.5. We see that this distribution is enhanced around m-^n = 0.4 GeV due to 
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Figure 3: The ratio ^^ — /^^ — in two regions of 7X1^^.,^, around 0.4 GeV (left panel) and 
around 1.3 GeV (right panel). 



the S-wave contribution and around m^. 



1.3 GeV due to D-waves. To illustrate the 



dN' 



-I 



dN' 



at different values of 



enhancements in these regions we plot in Fig. |^ the ratio ^^ / ^^ 
Bjorken x. From these figures we see that the most favorable regions to observe isoscalar 
channel in vr+vr" productions (and hence the gluonic component of 27rDA) are xbj > 0.3 
and mTTTT around 0.4 GeV (S-wave) and 1.3 GeV (D-wave). 



On Fig. we show the estimates of the ratio of — 



do-" 



dt It— tmin 

1.4 GeV 



to Tt \t-t 

dt K — tmin 

We see that in the valence 



integrated 



over rrinn from the threshold m^7r = 2mT, to m^jr 

region xbj = 0.1-0.4 the cross section for tt^tt^ production is a few percent compared to 

that for TT+TT^ production (which is usually called p production). 

Let us note that the most promising way to determine the 1 = production amplitude 
is through its interference with the 1=1 amplitude. Experimentally this can be done 
studying the angular distributions of the produced 7r"'"7r~ pairs. For example, the intensity 
density (cos 6) {6 is the polar angle of the pion momentum in the CM frame with respect 
to the direction of the pion pair momentum q') directly probes the interference of the 
1 = and the 1 = 1 channel. It would be zero if the tt+tt^ pairs were produced purely 
in the 1 = 1 state. So the observation of a nonzero (cos 6) would clearly indicate the 
production of pion pairs in the isoscalar state. The corresponding expressions can be 
obtained easily using Eqs. (p!9| , ^0D . The full expression is rather tedious, here we give the 
simplified expression for the intensity density (cos 6) at 777,7^,7 close to the mass of the p 
meson: 



(cos0) 



1 



3|F^(mp)| Nf + ACp 



Im 



21: 



2/+ + It 



^foirrip) - 2/2(771,; 



(46) 



Numerically we get for {cos 9) values of about -4%, -7%, and -10% at 777^^ = nip and 
Xbj = 0.1,0.2, and 0.3 respectively. At lower invariant masses (cos^) is larger. At 7777r7r = 
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Figure 4: T/ie raizo of -^^ — |t=tmin ^^ ~^~dt — l*=*min integrated over m.,^T^ from the threshold 
to 1.4 GeV as a function of Bjorken x. 

0.5 GeV for example the values of {cos 9) are -18%, -27% and -33% for xbj = 0.1, 0.2, and 
0.3 respectively. Detailed estimates of various intensity densities will be done elsewhere. 

Summary 

We have calculated hard exclusive two-pion production at the leading twist level to leading 
order in a^. The amplitude is expressed in terms of skewed parton distributions (SPD) 
and two-pion distribution amplitudes (27rDA). For the 27rDAs we used the asymptotic 
expressions. The comparison of the cross sections for the production of isospin one and 
isospin zero pion pairs shows that for small Bjorken x {xbj < 0.1) the isospin-one channel 
dominates, whereas in the valence region [xbj > 0.3) at certain invariant two-pion masses 
the contribution of isospin zero production is also sizeable. We showed that the isospin 
zero amplitude is dominated by the gluon 27rDA. Hence the cross section for isospin zero 
pion pairs gives access to the gluon distribution in the pions. To extract the isospin 
zero channel measurements in the quark valence region are needed. This extraction is 
facilitated if vr^vr^ pairs can be detected since the isospin one channel contributes only to 
TT+TT" pairs. 

Let us note that the general expression for the leading twist amphtude Eq. @ can 
be used to study other channels like KK, ttK, etc. By crossing symmetries this allows 
to obtain information on the partonic distributions of the produced particles. This is 
particularly important when the particles are not available as targets in deep inelastic 
reactions. The exclusive two-particle production can be accessible for measurements e.g. 
in the COMPASS experiment ^. 

We are aware of the fact that higher twist corrections could be large. They have not 
been analyzed yet but there is some hope that they cancel at least partially in the ratios 
of cross sections we discussed. 
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